function Mat = vdm11(d)
I=d:-1:0;
J=0:d;
IM = diag(I)*ones(d+1,d+1);
JM = diag(J)*ones(d+1,d+1);
Mat = (IM/d).^(IM').*(JM/d).^(JM');
IF = gamma(I+1);
JF = gamma(J+1);
A = factorial(d)*ones(d+1,d+1)*diag(1./(IF.*JF));
Mat = A.*Mat;